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Abstract —Massive Multiple-Input Multiple-Output (MIMO) 
is foreseen to be one of the main technology components in 
next generation cellular communications (5G). In this paper, 
fundamental limits on the performance of downlink massive 
MIMO systems are Investigated by means of simulations and 
analytical analysis. Signal-to-nolse-and-lnterference ratio (SINK) 
and sum rate for a single-cell scenario multi-user MIMO are 
analyzed for different array sizes, channel models, and precoding 
schemes. The impact of hardware impairments on performance 
is also Investigated. Simple approximations are derived that show 
explicitly how the number of antennas, nnmber of served users, 
transmit power, and magnitnde of hardware impairments affect 
performance. 

I. Introduction 

Massive Multiple-Input Multiple-Output (MIMO) systems 
are foreseen to be one of the main technology components 
in next generation cellular communication systems (5G). The 
basic idea with massive MIMO is to use a large number 
of antenna elements at the base station (BS)—an order of 
magnitude larger than used in current systems and much 
larger than the number of concurrently served user equipment 
(UE)—in order to achieve high spatial resolution and array 
gain. The high spatial resolution enables high system capacity 
by spatial multiplexing of UEs. The array gain enables high 
energy efficiency. Eurthermore, the averaging effects obtained 
by using a large number of antenna elements make the channel 
behave almost deterministically, which has the potential to 
simplify radio resource management. 

In order for massive MIMO to be economically feasible, the 
cost per antenna element and its associated radio and base band 
branches must be signihcantly less than in current systems. 
Therefore, the requirements put on such components must be 
less stringent than in current systems. Hence, massive MIMO 
can be seen as a paradigm shift from using a few expensive 
antennas to many cheap antennas JT]. 

Contributions: In this paper, a performance analysis 
of downlink massive MIMO is presented. The focus is on 
the impact of basic antenna and channel model parameters 
on system performance. Downlink signal-to-interference ratio 
(SINR) and sum rate for the matched hlte|^ (ME) and zero 
forcing (ZF) precoders are computed for a single-cell sce¬ 
nario using line-of-sight (LoS), independent and identically 
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^Sometimes also referred to as conjugate beamforming or maximum ratio 
transmission. 


distributed (IID) Rayleigh, and statistical ray-based channel 
models. Simulation results are presented together with simple 
analytical expressions where the throughput dependence on 
the system parameters is transparent. Furthermore, the impact 
of phase and amplitude errors in the precoding weights is 
analyzed. 

Previous works: A lower bound on the massive MIMO 
downlink performance has been derived in Q, under the 
assumption that the propagation channel is learnt through 
the transmission of orthogonal pilot sequences in the uplink. 
Hardware impairments in massive MIMO have been the sub¬ 
ject of several recent theoretical investigations. Most of these 
works model the hardware impairments as power-dependent 
Gaussian additive noise, which simplihes the throughput anal¬ 
yses 1^. More realistic multiplicative hardware models have 
been adopted in e.g., 0. Comparisons between the downlink 
throughput corresponding to measured propagation channels 
and ideal channel models have been reported in e.g., Q. 

II. System Model 

We consider a single cell with no external interference. The 
downlink signals received by K co-scheduled UEs, each with 
a single antenna, served by a BS with M antenna elements 
are modeled by the following equation 

y = vTHWx-l-e. (1) 

Here, y = [ t/i • • • yK is a AT x 1 vector containing the 
signals received by the UEs, P is the BS transmit (Tx) power 
and H = [ hf ... is the K X M channel matrix, 

with hfe being the 1 x M channel vector from the M BS 
antenna elements to UE k. Furthermore, W = [wi... Wif] 
denotes the M x K precoding matrix, x the K x 1 vector 
of transmitted symbols, and e a AT x 1 zero-mean complex 
Gaussian vector with covariance matrix NqI modeling the 
receiver noise in the UEs. 

When increasing the number of antenna elements in the 
analysis, the Tx power will be scaled as 1/M in order to 
compensate for the increased array gain. As a result, a constant 
signal-to-noise ratio (SNR) operating point is maintained. 
More precisely, the Tx power for a given number of antenna 
elements is set so that a certain target SNR, SNRt, is obtained 
in the interference-free (IF) case, i.e., 

P = NoSNRt/M. (2) 


The BS antenna array is assumed to be a uniform lin¬ 
ear array (ULA) in the horizontal plane with 0.6A element 
separation, where A is the wavelength. The radiation pattern 
of a single element is modeled according to Q with 90° 
azimuth half-power beam width. Mutual coupling between 
array elements is ignored. The UE is assumed to have a single 
isotropic antenna. Polarization is not modeled. 


III. Perfect Channel State Information 

In this section, the system throughput is analyzed assuming 
that the BS has perfect channel state information (CSI). The 
SINR and sum rate are calculated for ME and ZE precoding 
schemes. Some analytical results are first derived for an IID 
Rayleigh channel model and then compared with other channel 
models by means of simulations. 

The SINR for UE k for a given channel realization H is 
given by 

+ iVo' 

The maximum achievable rate for UE k is then given by 

= E [log 2 (1 -f 7fc)] (4) 


where the expectation is taken over the channel realizations. 
The sum rate is obtained by summing the rates of all concur¬ 
rently served UEs according to 

R = Y.Rk- ( 5 ) 

k 

A. IID Rayleigh 

In this section, simple approximations of the rate in Q are 
derived under the assumption of an IID Rayleigh channel, i.e., 
a channel for which H has IID zero-mean complex Gaussian 
elements. The expectation in Q is difficult to compute ana¬ 
lytically. However, as will be shown by simulations later in 
this paper, the following approximation is accurate 


i?fc«log2(l + SINRfe), (6) 


where 


SINRfc = 


PE[|hfcWfcp] 

PYd=i E[|hfeW2|2] -h Nq 


( 7 ) 


Analytical approximations of the SINR for ME and ZE pre¬ 
coders are derived below. These can be used together with 
(|^-(|^ to estimate the sum rate in the cell. 

1) MF: The matched filter precoding vector for UE k is 


Wfc = hf/||hfc||, 


(8) 


where the normalization by ||hfc|| assures the same transmitted 
power in all channel realizations and equal power allocation 
to all UEs. The expectations in Q are then given by 


E[|hfcWfcp] = E 


Ihfchf / 


= E[||hfef]=M, (9) 


and 


Hence, 


E[|hfcw,f] =E 


SINRfe = 


Ihfchf/llh 


PM 


= 1 . 


P{K - 1) + iVo ■ 

Using the power normalization in (|^ we obtain 

SNRt 


SINRfc = 


1-f SNRt(A:- 1)/M’ 


( 10 ) 

( 11 ) 


( 12 ) 


2) ZF: With equal power allocation to all UEs, the ZE 
precoding vector is 


Wfc = (Hl)fc/y||(Ht)fe||^, (13) 

where (Hl)fc denotes the fc-th column of Hi, Ij-lji? the 
Erobenius norm, and Hi = H^(HH^)“^. In order to 
simplify the analytical calculations, the ZE precoding matrix 
is approximated by 

W = cH1, (14) 


where c is a normalization constant obtained by solving 

E[||W||^]=iT. (15) 


This normalization makes the average transmitted total power 
to all K UEs equal to K, but does not guarantee same power in 
all channel realizations and equal power allocation to all UEs. 
However, for an IID Rayleigh channel the difference turns out 
to be small when the number of antenna elements is large. 

Under the assumption of an IID Rayleigh channel, HH^ 
is a K X K central complex Wishart matrix with M degrees 
of freedom and covariance matrix equal to the identity matrix. 
It follows then that ||^ p. 26] 

E[Ti {{nn^)-^}] = K/{M - K) (16) 


and the normalization constant is obtained according to 


c^E ||Ht||t =K ^ c = yjM - K. 

M W r 

The received signal according to Q is then given by 
y = ^P(M- A:)HH^x -f e = ^/P{M - K)x F e. 


(17) 

( 18 ) 


Hence, 


SINRfc = P{M - K)/No = SNRt(l - K/M). (19) 

3) Simple rule of thumb: The SINR approximations in o 
and ( [T9l l can be used to derive a simple rule of thumb on 
how many antennas are needed to reach a certain performance 
for a given number of co-scheduled UEs. Eor example, the 
number of antennas needed to reach an SINR that is 3 dB away 
from the IE SINR, SNRt, can be obtained by setting SINR = 
SNRt/2. This leads to the following number of antennas for 
ME 


M = (AT-l)SNRt, (20) 


and for ZE 


M = 2K. 


( 21 ) 
















4) Comparison between analytical and simulation results: 
Figured shows a comparison of the analytical approximations 
in and ( [T9l l with the the exact rate expression in 

Q obtained by Monte Carlo simulationsj^ The plot shows 
the average sum rate vs. number of antenna elements for 
K = IQ and SNR^ = 10 dB. Clearly, there is excellent 
agreement between the analytical and simulation results, also 
for moderately sized arrays. These results provide empirical 
support to the approximations made in deriving the analytical 
expressions. 



Fig. 1. Average sum rate vs. number of antennas for an IID Rayleigh channel 
using MF and ZF precoders. Comparison between analytical and simulation 
results with K = 10 and SNRt = 10 dB. 

B. Comparison of Channel Models 

The IID Rayleigh model is reasonable when there is rich 
scattering around the BS and the UEs. However, in many cases 
the environment around the BS is such that there is spatial 
coiTelation among the BS antenna elements. We investigate 
the impact of spatial correlation on massive MIMO by using 
two types of correlated channel models; 

1) LoS channel model: We assume that there is a single, 
free-space planar wavefront from the BS to each UE. Eor a 
given azimuth angle of departure (AoD), the model is purely 
deterministic. However, the AoDs to different UEs are drawn 
from a uniform distribution over the interval [-60°, 60°]. 

2 ) Statistical ray-based channel models: Two different ray- 
based models are investigated: the 3GPP spatial channel model 
(SCM) 0 and the ITU urban macro (UMa) model Q with 
indoor UEs added as described in H). 

The IID Rayleigh and LoS channel models represent corner 
cases in terms of angular spread, with the IID model being 
spatially white and the LoS model having zero angular spread. 
The ray-based models lie in between with a mean angular 
spread of 15° in the SCM model and 14° and 26° mean 
angular spread for LoS and non-LoS UEs, respectively, in the 
ITU UMa model. 

The SCM and ITU channel models also include models for 
path loss. In order to isolate the impact of path loss and spatial 
correlation on massive MIMO performance and to be able to 
compare with the unit gain IID Rayleigh channel, the path loss 
is first removed from channel matrices generated by the SCM 
and ITU models. 

^In the Monte Carlo simulations of ZF performance, the precoding vector 
in (H) has been used. 


Figure shows a comparison of performance with the 
different channel models when the channel gain has been 
normalized to one. With MF, the best performance is achieved 
with a LoS channel, while the LoS channel gives the worst 
performance when using a ZF precoder. The MF behavior can 
be explained by the low sidelobes attainable by a large antenna 
array, which result in good interference suppression in a LoS 
channel. The cause of the ZF behavior is that there is a gain 
penalty when two UEs are separated by less than a beam width 
in a LoS channel. Placing a null in the direction of one UE 
will, in this case, give a large gain drop to the other UE. In an 
uncorrelated channel, however, it is unlikely that two channel 
vectors are almost parallel. 

The other channel models give very similar performance for 
both the ME and ZE precoders. It is interesting to note that 
the SCM and ITU channel models give almost identical per¬ 
formance as an IID Rayleigh channel. A possible explanation 
to this is that different UEs will obtain different realizations 
of the channel rays. The only correlation between UEs is 
via the correlation between the large-scale parameters, such 
as angular spread, delay spread, and shadow fading. User- 
common clusters are not captured by these models. This is 
something that may be important when evaluating multiuser 
MIMO performance. 




Fig. 2. Average sum rate vs. number of antenna elements for different channel 
models, top: MF, bottom: ZF. 


In order to investigate the impact of path loss difference 
between UEs on performance, Eigure shows the average 
sum rate relative the corresponding IE case as a function 
of the number of antenna elements using the ITU UMa 
model with and without path loss. The results show that ME 
relative performance is reduced significantly when using path 
loss in the channel model, while the impact on ZF relative 
performance is weak. 































Fig. 3. Average sum rate relative IF vs. number of antennas for different 
precoders with the ITU UMa channel with and without path loss 


IV. Imperfect Channel State Information 

Perfect CSI at the BS cannot be achieved in reality so the 
results in the previous section should be interpreted as upper 
bounds on performance. In this section, we analyze the impact 
of imperfect CSI using a simple model of hardware impair¬ 
ments. We model imperfect CSI as a phase and amplitude error 
applied to the true channel. More specifically, the channel used 
in the downlink transmission is modeled as 

hkm “ (1 “f — ^mhkm (22) 

where hkm is the perturbed downlink channel coefficient 
between user k and antenna m, hkm is the corresponding true 
channel coefficient, and — (1 + am)eP^"'■ The amplitude 
error, am, and phase error, (j>m, are assumed to be independent, 
zero-mean Gaussian random variables with variances and 
cr^, respectively. In the remainder of the paper, we shall refer 
to a a as amplitude error, and measure it in dB, and to as 
phase error, and measure it in degrees. 

The adopted hardware impairment model is not a model of a 
particular component. Rather, it captures the aggregated effect 
of all errors in the system. The polar form of this aggregated 
impairment model is supported by the fact that two of the 
largest impairments, i.e., the power amplifier distorsion GO) 
and the oscillator phase noise G3 are multiplicative. It has 
been shown in Q that system performance predictions based 
on this model are in good agreement with the ones based on 
more refined hardware models. 

1) MF: In the case of errors, the MF precoding vector is 
given by the perturbed channel vector 

Wfc = hf/||hfe||. (23) 

To simplify calculations, we approximate the norm by its 
expected value and instead use the following precoding vector 

Wfe = hf/-\/E[||hfc||2]. (24) 

As we shall later show, this is an accurate approximation for 
large M. Using ( |22l l we obtain 

M 

E[||hfc||2] = E[^ \{l+am)e^^-hkm\^] = (25) 

m—1 





■ M 

E 

IhfchfP 

= E 

1 'y ] hkmhkm^rn] 

_ m—1 


= E 


M MM 


TTl—l 


m—1 '^=1 


Since E[|z|‘^] = 2cr'* for a complex Gaussian random variable 
z with variance a^, we have that E[|him|"^] = 2. To compute 
E[em], we use that the characteristic function, dehned as 
i/)(s) = E[exp(jsX)], of a Gaussian random variable X with 
mean fj, and variance is given by tp{s) = « /2 gy 

setting s = 1, we get 

(26) 


Since, by assumption, is a zero-mean Gaussian random 
variable with variance cr^, we obtain 

E[e„] = E[1 + am] . (27) 

The expectation in ( [26] l then simplifies to 

E[|hfchf n = 2M(1 -f a^) + M{M - (28) 

which, for large M, can be approximated by 

E[|hfchf|2] (29) 


Finally, 


E[|h,.h"|^] =E 


M 

E 

m—1 




= M{l + al). (30) 


The SINR approximation in (j?]) is then given by 

SNRt 


SINRfc = 


(31) 


l-fCT2i + sNRt(R:-l)/M' 

Hence, the error-free SINR in 0 is reduced by a factor 
exp(—CT^)/(1 + in the presence of phase and amplitude 

errors. The factor exp(—cr^) reflects that a phase error causes 
the ideal and perturbed channel vectors not to be parallel. An 
amplitude error does not destroy the alignment, but yields 
a gain reduction by a factor 1/(1 + a^) due to the weight 
normalization in ( |24l i which is needed to assure conservation 
of energy. If the errors are small, the error factor in ( |3l] l can 
be approximated by Taylor expansions according to 


l + o-? 


1 + cr^ 


1 + ^/ 


(32) 


where cr^ = + cr^ is the total error variance. Hence, for 

small errors, the SINR degradation depends only on the sum 
of the variances of the phase and amplitude errors. The factor 
in ( [32l l is the same as the gain reduction caused by phase and 
amplitude errors in phased arrays G3- Note that this SINR 
reduction does not depend on the number of antenna elements. 
Hence, it remains even in the limit M ^ oo. However, with 
the Tx power normalization used in this paper, the system will 
asymptotically be noise limited; so the asymptotic SINR loss 
can be compensated for by an increase in Tx power. 


and 

























2 ) Simple rule of thumb: A rale of thumb for the required 

number of antennas to be 3 dB from the IF SINK is easily 
derived also in the case of phase and amplitude errors. Using 
the approximations in ( |3T] i and ( |32l i, the following expression 
is obtained ^ 

M= ^i^(A:-l)SNRt. (33) 

Hence, the required number of antennas is increased by a 
factor (l+cr^)/(l—(T^) compared to the error-free case in ( |20| . 

3) Simulation results: Figure shows a comparison be¬ 
tween simulation results using an IID Rayleigh channel model 
and the approximations pT])-(|3^ when the standard deviation 
of the phase and amplitude errors is 20° and 1 dB, respectively. 
The agreement between simulation and analytical results is 
excellent. Since analytical expressions for imperfect CSI have 



Fig. 4. Average sum rate vs. number of antennas for an IID Rayleigh channel 
using an MF precoder with 1 dB amplitude and 20° phase errors . Comparison 
between analytical and simulation results. K = 10 and SNRt = 10 dB. 

been derived only for the MF precoder some simulation results 
including also the ZF precoder are now presented. Since the 
IID, SCM, and ITU UMa channel models give very similar 
results only results for the IID model are given here. 

Figure shows average sum rate relative the error-free case 
vs. amplitude error when there is no phase error and vice 
versa for MF and ZF precoders with 20, 100, and 500 antenna 
elements, respectively. The results show that the ZF precoder 
is more sensitive to errors than MF. It can also be seen that 
the impact of errors on the relative sum rate decreases as the 
number of antenna elements is increased. 

V. Conclusions 

In this paper, a performance analysis of single-cell, down¬ 
link massive MIMO has been presented. The results show that 
the 3GPP SCM and ITU UMa channel models give similar 
performance predictions as an IID Rayleigh channel if the path 
loss is the same for all UEs. Including path loss differences 
between UEs gives a significant performance reduction for 
the ME precoder, whereas the ZE precoder only gets a minor 
degradation. 

Eurthermore, the impact of phase and amplitude errors 
on performance was analyzed. Analytical analysis, based on 
approximations which are validated by simulations, show that 
these errors give an SINR loss that is independent of the 
number of antenna elements for the ME precoder in an IID 
Rayleigh channel. However, with the Tx power normalization 
used in this paper, the system will asymptotically be noise 


MF ZF 



MF ZF 



Fig. 5. Average sum rate relative the error-free case vs. amplitude (top) and 
phase (bottom) error for MF (left) and ZF (right) precoders. A" = 10 and 
SNRt = 10 dB. 

limited. Therefore, for large arrays, this loss can be com¬ 
pensated for by an increase in Tx power. It was shown by 
simulations that the ZE precoder is more sensitive to phase and 
amplitude errors than the ME precoder. Some simple analytical 
approximations for SINR and sum rate have been derived for 
the IID Rayleigh channel model. These approximations show 
how different parameters impact system performance. 
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